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(1  eV  to  20  eV)  electrons  in  silicon  dioxide  with  and  without  electric  fields; 

(2)  comparison  testing  of  one-dimensional  discrete  ordinates  transport  codes  for 
charged  and  neutral  particles  with  analytical  benchmark  calculations;  (3)  Monte 
Carlo  transport  calculations  for  low  energy  electrons  in  the  presence  of  strong 
electric  fields  in  silicon  dioxide;  (4)  Monte  Carlo  transport  calculations  for  high 
energy  (MeV)  electrons  in  insulating  materials  with  Internal  electric  fields. 
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I.  Introduction 


During  the  period  from  September  1,  1989  through  July  22,  1992,  the  execution  of 
contract  No.  F19628'89-C-0189  was  directed  toward  the  performance  of  mathematical 
analysis,  research  and  computer  program  development  for  the  simulation  and  analysis  of 
charged  particle  transport  in  irradiated  and  non-irradiated  solids. 

This  report  is  organized  into  five  sections  in  addition  to  this,  the  introduction. 
Section  II  is  a  discussion  and  presentation  of  research  in  discrete  ordinates  transport 
calculations  of  low  energy  (  ~  1-20  eV)  electrons  in  silicon  dioxide  with  and  without 
applied  electric  fields  such  as  are  present  in  microelectronic  devices.  Section  III  is  a 
brief  description  of  research  activity  in  the  development  and  application  of  analytical 
benchmark  calculations  for  one-dimensional  transport  in  the  evaluation  of  the  accuracy 
of  approximate  discrete  ordinates  methods.  Section  IV  is  an  account  of  our  Monte 
Carlo  transport  calculations  for  low  energy  electrons  in  strong  electric  fields  in  silicon 
dioxide.  These  calculations  were  developed  for  the  purpose  of  describing  and  predicting 
the  extent  and  significance  of  the  impact  ionization  process  in  microelectronic  device 
materials.  Section  V  presents  the  results  of  another  set  of  Monte  Carlo  electron 
transport  calculations  which  we  developed  for  describing  the  transport  of  MeV  electrons 
in  insulating  materials  with  internal,  polarization-induced  electric  fields.  Section  VI, 
the  final  section,  is  the  reference  list. 
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n.  Discrete  Ordinates  Electron  Transport  Calculations 

1.  Introduction 

In  earlier  research  efforts  involving  the  method  of  discrete  ordinates  (S;^)  for  the 
solution  of  electron  transport  problems,  we  had  made  extensive  use  of  the  transport 
code  ONETRAN^^%  as  reported  in  Ref.  2.  As  the  research  activity  continued,  however, 
it  became  clear  that  creation  of  a  new  computer  program  with  the  ONETRAN  linear 
discontinuous  and  diamond  differencing  algorithms  would  be  highly  useful.  It  was  our 
intention  to  apply  the  S;v  method  to  a  variety  of  electron  transport  problems  covering  a 
wide  range  of  electron  energies  and  scattering  processes,  from  the  eV  range  in  which 
electron-phonon  scattering  is  important  to  the  MeV  range  in  which  the  Mott  scattering 
model  is  applicable.  It  became  apparent  to  us  that  if  the  Sjy  method  were  to  be 
effectively  applied  to  electron-phonon  scattering  in  microelectronic  device  materials,  it 
would  be  necessary  to  build  into  the  S^v  program  a  capability  for  inclusion  of  applied 
electric  field  effects,  since  the  voltages  typically  applied  to  silicon  devices  produce  large 
internal  electric  fields.  Algorithm  modifications  such  as  this  would  be  difficult,  if  not 
impossible,  to  accomplish  with  the  ONETRAN  code  without  encountering  serious 
problems  that  could  arise  from  programming  error  and/or  algorithm  incompatibility. 

The  applicability  of  SNCODE  is  restricted  to  slab  geometry.  The  sacrifice  in 
generality  has  resulted  in  a  program  that  is  simpler  than  ONETRAN,  which  treats 
cylindrical  and  spherical  geometries  as  well.  Another  departure  from  the  ONETRAN 
method  is  that  SNCODE  computes,  preserves  and  uses  the  angular  components  of  the 
flux,  while  ONETRAN  does  the  same  with  the  Legendre  series  expansion  coefficients  of 
the  angular  flux.  The  actual  angular  components  are  computed  for  the  spatial 
differencing  procedure  and  then  discarded.  The  ONETRAN  method  is  more 
computationally  efficient  when  the  scattering  cross  sections  are  given  in  terms  of  their 
Legendre  expemsion  coefficients,  which  is  usually  the  case.  While  the  two  methods  are 
equivalent,  it  is  not  convenient  to  apply  the  ONETRAN  approach  when  dealing  with 
applied  electric  fields  oriented  along  a  single  direction.  In  the  following  sections  we  will 
discuss  various  applications  of  SNCODE  to  :  l)electron-phonon  scattering  with  and 
without  applied  electric  fields;  and  2)benchmark  calculations  and  verification  of 
differencing  algorithm  accuracy. 

2.  Low-Energy  Electron  Transport  Calculations 

2.1  High-Resolution  (Energy)  Calculations  of  Electron-Phonon  Scatter  in  Silicon  Dioxide 

The  discrete  ordinates  transport  method  was  used  to  calculate  electron  energy 
spectra  resulting  from  electron-phonon  scatter  in  SiOj  under  zero  electric  field 
conditions.  SNCODE  was  installed  and  run  on  the  Phillips  Laboratory  CRAY-2.  The 
availability  of  the  supercomputer  provided  an  opportunity  to  make  very  high  resolution 
calculations  of  these  electron  energy  spectra.  A  cross  section  set  Wcis  generated  for  2000 
energy  groups  of  uniform  width  over  the  interval  0-10  eV  (AE=0.005  eV).  This  allowed 
full  recognition  of  the  discrete  optical  mode  energy  transitions  (.068,  .153  eV)  and  a 
reasonably  accurate  representation  of  the  smaller  acoustic  mode  energy  changes 
(maximum  energy  transfer  ~  0.02  eV).  A  5  eV  electron  source  was  positioned  at  the 
center  plane  of  a  Si02  slab  of  thickness  100  nm.  The  source  was  directed  to  the  right 
side  of  the  center  plane  with  a  cosine- squared  (^^;  see  insets  of  Figs.  1  and  2)  angular 
distribution.  The  2000  group  run  was  made  using  an  energy  group  in-scatter  source 
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error  toleramce  of  t  =  10“^.  That  is,  a  number  of  outer  (energy)  iterations  were  made 
until  the  group  with  largest  error  in  in-scatter  source  reached  a  value  <  (g.  For  this 
case,  the  required  number  of  outer  iterations  was  29.  The  highest  energy  group 
containing  up-scattered  electrons  (in-scatter  source  given  the  above 

convergence  condition,  was  group  883,  corresponding  to  ~  5.6  eV.  The  resulting  energy 
spectra  show  the  optical  mode  transition  lines  at  the  low  energy  end  superimposed  on 


Fig.  1.  Electron  energy  spectrum  at  +40  nm  in  SiOj  under  zero 
field  conditions.  Electron  source,  5  eV,  located  at  center  plane  of 
lOOnm  slab  is  distributed,  directed  to  the  right  of  center  plane. 
Spectrum  structure  results  from  electron-phonon  interactions. 
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Energy  (eV) 


Fig.  2.  Electron  energy  spectrum  at  -40  nm  in  Si02  under  zero 
field  conditions.  Electron  source,  5  eV,  located  at  center  plane  of 
lOOnm  slab  is  distributed,  directed  to  the  right  of  center  plane. 
Spectrum  structure  results  from  electron- phonon  interactions. 


the  largely  acoustic  mode  generated  smooth  spectrum.  As  evidenced  by  the  spectra 
shown  in  Figs.  1  and  2,  the  Syy  method  provided  a  degree  of  energy  spectrum  resolution 
not  obtainable  with  the  Monte  Carlo  method.  The  two  spectra,  shown  at  +  40  nm  (Fig. 
1)  and  -40  nm  (Fig.  2),  are  similar.  The  amount  of  electron  scattering  at  these 
distances  is  sufficient  to  nearly  eradicate  any  differences  due  to  source  directionality. 
The  spacing  of  the  peaks  is  consistent  with  the  values  of  the  optical  mode  energy 
transitions  given  above. 
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2.2  Discrete  Ordinates  Treatment  of  Low  Energy  El*  'tron  Trttnspor-  in  the  Presence  of 
an  Electric  Field 

Modifications  were  mad-  o  the  S;v  program,  SNCODE,  to  account  for  the  effects 
of  the  presence  of  an  electric  field  on  electron  transport  in  solids.  The  objective  of  this 
effort  was  to  model  low  energy  electron  transport,  electron- phonon  scattering,  and 
production  of  secondary  electrons  due  to  impact  ionization,  in  the  oxide  layer  of  a 
silicon  device  under  typical  operating  conditions,  t.e.  applied  voltage  «  5  V.  Much  of 
our  method  for  the  treatment  of  electric  field  effects  was  based  on  the  earlier 
formulation  of  Wienke^^^  in  which  he  incorporated  a  Lorentz  force  term  into  the 
Boltzmann  equation  (no-relativistic)  for  electrons. 

The  multigroup  non-relativistic  Boltzmann  transport  equation  (in  one-dimension) 
for  group  g  electrons  in  the  presence  of  an  electric  field  is 

Qfg  "I"  Qag  (^) 

where 

4>g(x,ii)  is  the  group  g  flux, 

X  and  /i  are  the  position  and  velocity  direction  cosine,  defined  in  the  usual  way, 

<7^  is  the  total  cross  section, 

‘i  is  the  electric  field  strength, 
e  is  the  electronic  charge, 
ttiq  is  the  electron  rest  mass, 

Vg  is  the  velocity  of  group  g  electrons, 

Q fg  is  the  fixed  source  for  group  g, 
and  Qgg  is  the  in-scatter  source  to  group  g. 

Our  approach  differed  from  that  of  Wienke  in  that  we  began  with  the  Boltzmann 
equation  rather  than  the  Fokker-Planck  equation.  SNCODE  is  a  Boltzmann  solver  and, 
as  will  be  demonstrated  elsewhere  in  this  document,  operates  correctly  for  applications 
in  which  no  electric  field  is  present.  We  treated  the  velocity  gradient  electric  field  term 
in  the  same  manner  as  Wienke. 

Since  our  multigroup  discrete  ordinates  equation  solves  for  the  flux  as  a  function 
of  energy,  i.e.  <t>g{x,n)  as  ^(x,n,Eg),  it  was  necessary  to  re^express  the  velocity  gradient 
term  in  terms  of  the  energy  variable.  The  electric  field  8  was  assumed  to  be  directed 
along  ~r.  With  the  electron  velocity  given  by  "v,  and  the  velocity  direction  cosine  with 
respect  to  the  field  direction  given  by  /i  (see  Fig.  3), 


Fig.  3.  Vector  diagram  showing  relative 
orientations  of  electric  field  and  electron  motion. 
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The  following  useful  identity  for  the  velocity  gradient  operator  is  given  in  Ref. 


3  as 


Defining 


then 


(3) 

(4) 


jj^6  •  =  1.  +  ^  . 

(5) 

The  procedure  for  incorporation  of  the  velocity  derivative  term  L  into  the  multigroup 
Boltzmann  will  be  outlined  first,  and  that  for  the  angular  derivative  term  will  be 

given  second. 

The  operator  L  can  be  simplified  as  follows: 

•  ^  _  e  y  V  d4> 

E4>  -  rh(,t;G  •  V  Qy 

_  e  »~r  'v  dtj) 

-  m^T-  V  ^ 

_  c/i 

“  rnyv  ^  • 

(6) 

Our  multigroup-S^v  calculation  requires  that  the  flux  be  expressed  in  terms  of  energy 
rather  than  velocity.  For  non-relativistic  calculations,  the  velocity  derivative  can  be 
directly  expressed  in  terms  of  energy  derivatives. 

d  dE  d  p  1  2 

(7) 

thus 

(8) 

O  1 

and  L4t  -  eSfi^  . 

(9) 

Eqn.  (1)  can  then  be  restated  as 

dib  Q 

=  Q/ff  +  • 

(10) 

In  finite  difference  terms,  the  energy  derivative  term  becomes 


I 

I 


e8/i  -  for  ^  >  0, 

eSfi  for  ^  <  0  . 


(lla) 

(llb) 


The  /i>0  case  corresponds  to  energy  upscatter  by  the  field,  and  the  /i<0  case 
corresponds  to  energy  downscatter.  In  the  former,  the  positive  velocity  direction  cosine 
indicates  that  the  electron  is  moving  with  the  field  direction,  while  in  the  latter,  the 
•ctron  motion  is  in  opposition  to  the  field  direction.  If  we  order  the  energy  groups 
such  that  E  -y  >  Eg  >  Eg_^,^,  and  designate  the  energy  group  width  at  Eg  cis  A£'g, 
then  the  field  acceleration  terms  become 
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for  /I  >  0, 


(11a) 


L  ss  e8/i  ^  — 

./ 

I  m  -  e6/i  —  for  p  <  0 


(11b) 


The  finite  difference  (in  the  energy  variable)  form  of  Eq.  10  then  becomes 

^  +  i  +  Q/ff  +  Qtg  (;^  >  0)  .  (12a) 

=  "1^  +  Q/tf  +  Q*a  (/^  <  0)  •  (121>) 

a  a 

In  Eqs.  12a, b  the  terms  involving  4>g^i  have  been  transposed  to  the  right-hand  (source) 
side  and  can  be  regarded  as  positive,  angularly  dependent  contributions  to  the  in-scatter 
source.  The  e^nfAEg  term  on  the  left  side  has  the  effect  of  adding  a  positive  angularly- 
dependent  term  to  the  total  cross  section. 


The  angular  derivative  term. 


can  be  written 


_  e 


iHqV  c/1 


which  has  the  same  form  as  the  angular  part  of  the  Fokker-Planck  r..*a,tor  and 
represents  a  within-gioup  angular  redistribution  (trajectory  bending)  by  the  electric 
field.  Because  electron-phonon  interactions  tire  largely  characterized  by  wide  angle 
scatters,  the  extended  P12  transport  correction  adequately  handles  the  scattering 
angular  distribution.  After  a  small  number  of  scatters,  the  flux  is  not  highly 
anisotropic,  and  the  angle  derivative  can  be  evaluated  with  a  flux  Legendre  expansion  of 
order  12;  assuming  that  12 

«  £  (<  +  5)  (‘i  P,M  (H) 

1  =  0 

then  making  use  of  the  identity 


yields 


where  the  flux  moments  <l>*  are  obtained  with  well-known  expression 


+ 1 

4'g  =  j  d/4  4>g(fi)Pl{fi)  . 


With  the  above  expression  for  Jt,  the  Boltzmann  equation  is  then  given  by 

+[4  ^  ^  ^  ijiUa  =  =*=  ll;  ^±1  +  Q/tf  +  (/^  >  0)  ,  (16) 
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where  the  ±  denotes  +  for  /i  >  0  and  -  for  p  <  0  . 

In  the  above,  it  is  assumed  that  energy  transfers  due  to  electric  fields  occur 
between  adjacent  groups,  the  continuous  approximation  (CDSA).  Eq.  16,  the 
multigroup  transport  equation,  is  solved  in  the  usual  manner  by  the  method  of  discrete 
ordinates.  The  spatial  differencing  scheme  used  is  linear- discontinuous,  and  the  discrete 
ordinate  angles  are  chosen  as  Gauss  quadrature  ordina,tes. 


The  formulation  given  above  was  implemented  in  SNCODE,  and  the  results  were 
reported  in  Bef.  4.  In  this  work,  the  cross  section  parameters  for  the  multigroup 
ra'culation6  vcre  derived  from  the  in'^erse  mean  free  path  and  stopping  power 
expressions  of  \3hley^®’®l  The  width  structure  of  the  energy  groups  was  chosen  in  order 
to  malce  ose  of  averaged  quantities  (cross  section,  stopping  power).  As  illustrated  in  the 
graph  of  Fig.  4,  the  stopping  power,  averaged  over  all  individual  electron-phonon 
interaction  modes,  is  negative.  Thus  in  the  term  of  Eq.  16  above,  the  group 
transfer  cross  sections  for  electron-phonon  interactions,  in  Eq.  17  below,  were 

characterized  only  by  downscatcer. 


= 


^7-  i~*a 


\ds/^  mode* 

AE„ 


(17) 


Up-scatter  can  only  occur  via  electric  field  action. 


Fig.  4.  Individual  mode  and  averaged  electron  stopping  powers  for 
electron-phonon  interactions  in  SiOj. 
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2.2.1  Electron  Multiplication 

To  account  for  the  secondary  electron  production  which  can  occur  through  impact 
ionization  in  Si02  when  the  primary  electron  energy  exceeds  8.9eV,  we  utilized  as  set  of 
ionization  cross  sections  can  calculated  from  the  theory,  developed  by  Ritchie^'^  and 
Ashley'*l  In  our  calculations,  impact  ionization  was  treated  as  (e,  2e')  scattering. 
The  energies  of  the  primary  and  secondary  electrons  after  the  ionizing  collision  were 
assumed  equal  ?ind  were  given  as  follows: 

E,  =  i  (Ei  -  E„,),  (18) 

where  E,  =  initiail  energy  of  primary  electron, 
and  E/  =  final  energy  of  primary  electron, 

=  energy  of  secondary  electron. 

The  factor  of  1/4,  rather  than  1/2  arose  from  the  allocation  of  an  equtil  cimount  of  recoil 
energy  to  holes.  Figure  5  is  a  plot  of  the  ionization  cross  section,  shown  together  for 
comparison  purposes,  with  the  phonon  scatter  cross  section.  It  is  seen  that  at  lOeV,  the 
ionization  component  is  lower  thain  the  phonon  scatter  by  more  them  two  orders  of 
magnitude.  However,  in  the  12-14eV  range,  this  factor  is  neirrowed  to  approximately 
15.  Thus,  at  sufficiently  high  energies,  significant  impact  ionization  is  possible. 


Energy  (eV) 


Fig.  5.  Comparison  plot  of  electron-phonon  interaction  and 
impact  ionization  cross  sections  in  SiOj. 
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2.3  Results 

2.3.1  5eV  Electron  Source;  Electric  Field  Strengths  =  0  and  10®V/cm 

S4  calculations  were  made  with  SNCODE  for  a  5eV  isotropic  source  of  electrons, 
uniformly  distributed  over  the  interval  40-60  nm  in  the  Si02  slab  (thickness  =  0.1  (xm) 
as  shown  in  Fig.  6.  A  uniform  electric  field  of  strength  8^  =  10®V/cm  was  assumed 
oriented  edong  the  z-direction.  The  energy  group  structure  consisted  of  100  uniformly 
spaced  groups  of  width  AE^  =  O.leV  spanning  the  energy  r«inge  O-lOeV.  The  spatial 
discretization  was  chosen  so  that  a  change  in  the  electron  position  of  one  Az  interval 
would  not  produce  a  field-induced  change  in  the  electron  energy  of  more  than  one  AE^  . 
This  is  an  essential  requirement  for  these  calculations,  since  it  is  assumed  that  energy 
transfers  due  to  the  electric  field  are  allowed  to  occur  only  between  adjacent  groups,  the 
“continuous  approximation” In  order  to  maintain  a  consistent  set  of  discretization 
parameters  in  these  Sff  calculations,  either  the  spatial  resolution  or  the  energy  group 
width  (or  both)  must  «Jso  be  increased  as  the  value  of  6^  is  increased. 

Fig.  6  shows  a  set  of  spatial  flux  profiles  at  8  energies  (1,  2,  3,  3.5,  4,  4.5,  5,  5.5 
eVl,  for  8,  =  10®V/cm.  The  flux  curves  shifted  in  the  direction  of  8^,  and  the  field- 
induced  up-scatter  is  evidenced  by  the  presence  of  flux  curves  at  energies  in  excess  of  5 
eV. 


Z  (micrometers) 


Fig.  6.  Electron  flux  profiles  at  10  energies  with  electric  field 
(6,=10*V/cm).  Unifornily  distributed  (.04-.06/im),  5  eV  isotropic 
electron  source  in  SiOj^^  . 
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The  effect  of  the  electric  field  on  the  distribution  of  electron  energies  is  also  shown 
in  Fig.  7  where  the  respective  energy  spectra  for  8^  =  0  and  8,  =  10®V/cm  are  plotted 
at  z  =  .09/im,  where  the  spectra  have  stabilized.  As  can  be  seen,  the  results  predict 
that  a  field  value  of  8^  =  10®V/cm  is  nearly  sufficient  to  promote  the  electron  energy 
above  the  8.9  eV  ionization  threshold. 


Fig.  7.  Electron  energy  spectra  for  uniformly  distributed  (.04- 
.06/iml  5  eV  isotropic  electron  source  in  SiOo  with  8,=0  and 
8,=10^V/cm.^^J 

2.3.2  8eV  Electron  Source;  Electric  Field  Strength  =  10  ®  V/cm. 

With  the  same  geometry  as  shown  in  the  Fig.  6  inset,  the  electron  source  energy 
was  increased  to  8eV,  a  value  slightly  below  the  ionization  threshold.  From  this 
starting  point,  electron  acceleration  due  to  the  10®V/cm.  field,  in  combination  with  the 
positive  contribution  of  the  phonon  absorption  terms  to  the  electron-phonon  interaction 
stopping  power,  were  sufficient  to  promote  the  source  electrons  to  energies  in  excess  of 
E™.  The  S;v  calculation  was  performed  using  140  uniformly  spaced  energy  groups  (AE 
=  O.leV),  100  uniform  spatial  steps  (Az  =  .001/im),  and  as  in  the  previous  case,  4 
discrete  ordinates.  With  the  outer  iteration  convergence  precision  of  10“^,  it  was  found 
that  15  outer  iterations  were  required  for  convergence. 
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Figures  8a  and  8b  show  the  electron  energy  flux  spectra  as  a  function  of  distance 
from  the  left  and  right  edges,  respectively,  of  the  source  band.  These  spectra  are  not 
absolute.  That  is,  at  each  vadue  of  z,  the  flux  profile  as  a  function  of  energy  is  shown. 
Thus,  while  the  low  energy  peak  (  ~  0.4eV)  in  Fig.  8a,  for  example,  appears  larger  than 
the  pefdc  at  ~  8eV,  the  opposite  is  actually  true  in  absolute  terms,  as  would  be  expected 
from  inspection  of  the  cross  section  curves  of  Fig.  5.  We  concluded  that  the  low  energy 
peak  in  Fig.  8a  was  primarily  due  to  downscatter  (phonon  emission)  and  deceleration 
caused  by  the  electric  field.  The  0.4eV  peak  occurs  at  a  distance  0.04  ;im  from  the 
source  band.  During  the  transport,  these  initially  back-directed  electrons  travel  in 
opposition  to  the  field  direction;  their  mean  free  path  increases  as  they  lose  energy;  and 
their  direction  of  travel  remains  “locked-in”  to  the  negative  direction.  The  electric  field 
alone  accounted  for  a  4eV  energy  loss  over  the  .04  ^m  distance.  On  the  right  side  of  the 
source  region  (Fig.  8b)  there  is  a  sharp  peak  at  the  8eV  source  energy,  with  a  high 
energy  tail  extending  to  ~  9.9eV.  Clearly  here,  energy  downscatter  still  dominates,  but 
the  energy  loss  rate  is  moderated  by  acceleration  due  to  the  field,  so  that  the  peak 
declines  gradually.  The  low  energy  peak  at  .04  fim  occurs  at  l.SeV.  This  may  be 
attributable  in  part  to  acceleration  by  the  field  of  both  downscattered  electrons  and 
secondary  electrons  arising  from  the  promotion  of  the  8eV  source  electrons  to  energies 
exceeding  E^„p. 

2.3.3  12eV  Electron  Source;  Electric  Field  Strength  =  2x10  ®  V/cm. 

Another  140  group  calculation  was  performed  for  a  12eV  source  with 
8  =2xlO®V/cm.  for  the  same  problem  geometry  as  above.  To  preserve  the  applicability 
of  the  “continuous  approximation”  with  regard  to  energy  gain/loss  due  to  the  field,  the 
resolution  of  the  spatial  discretization  mesh  was  doubled  (Ax=.0005/im)  beyond  that 
used  in  the  previous  case.  This  calculation  converged  with  fewer  outer  iterations  (12) 
than  were  required  for  the  10®V/cm.  case  (15).  Three  possible  reasons  for  this  are:  1) 
fewer  up-scatter  (phonon  absorption  events)  were  required  in  this  cailculation  to  reach 
the  14eV  upper  bound  due  to  the  higher  source  energy  in  this  calculation,  12eV  rather 
than  8eV;  2)  also  because  of  the  higher  source  energy,  secondary  electron  production 
began  with  the  first  outer  iteration;  and  3)  the  doubling  of  the  electric  field  force 
resulted  in  a  more  rapid  promotion  of  the  12eV  electrons  to  higher  energies. 

The  doubling  of  the  electric  field  strength,  in  combination  with  the  increased 
source  energy,  results  in  an  enhancement  of  low  energy  secondary  electron  production. 
Figure  9a  is  a  plot  of  the  electron  energy  spectra  vs.  distance  from  the  left  side  of  the 
source.  Unlike  the  10®V /cm.  case  shown  in  Fig.  8a,  there  is  a  small  peak  of  low  energy 
electrons  adjacent  to  the  source  region.  These  are  secondary  electrons  which,  when 
produced,  can  range  in  energy  from  0.75  to  1.5eV.  Similarly,  in  Fig.  9b,  the  spectra  on 
the  right  side  of  the  source,  a  low  energy  peak  near  the  source  region  appears,  also  in 
contrast  with  the  previous  case.  The  electric  field  strength  of  2xlO®V7cm.  was 
insufficient  to  overcome  the  phonon  emission  downscatter;  thus  the  secondary  electron 
energy  remained  low. 
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Flux  Profile 


Fig.  8a.  Electron  energy  flux  spectrum  profiles  vs.  distance, 
d(/im),  from  left  edge  of  source  region  in  Si02  (Fig.6  inset). 
Source  energy  =  8eV.  8,=10®V/cm.^^^ 


Fig.  8b.  Electron  energy  flux  spectrum  profiles  vs.  distance, 
d(/im),  from  right  edge  of  source  region  in  SiOj  {Fig.6  inset). 
Source  energy  =  8eV.  8  =10®V/cm.^^* 
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Fig.  9a.  Electron  energy  flux  spectrum  profiles  vs.  distance, 
d(/;<m),  from  left  edge  of  source  region  in  SiOj  (Fig.6  inset). 
Source  energy  =  12eV.  8,=:2xlO®V /cm.‘^* 


Fig.  9b.  Electron  energy  flux  spectrum  profiles  vs.  distance, 
d{fan),  from  right  edge  of  source  region  in  Si02  (Fig.6  inset). 
Source  energy  =  12eV.  6  =2xlO®V/cm/*^ 


2.4  The  “Continuous  Approximation”  Restriction 

In  the  calculations  described  in  Sections  2.2 -2.3,  we  found  that  imposition  of  the 
“continuous  approximation”  condition  restricts  the  applicability  of  the  Sj^  method  to 
electron-phonon  scattering  problems.  As  previously  stated,  with  this  condition  in  force, 
a  field-induced  change  in  energy  over  one  spatial  discretization  interval  must  not  exceed 
one  energy  group  width.  That  is, 

e&Ax  <  AEg  .  (19) 

For  field  strength  values  ~5  -  10  MV/cm  and  for  AEg  ~  0.1  eV,  the  range  of  Ax 
values  is  ~  1  -  2A.  T3'pical  SiOj  slab  thicknesses  for  these  calculations  were  in  the  500  - 
lOOOA  range.  A  low  energy  transport  calculation  for  0  -  20  eV  electrons  can  require 
~500  -  1000  armesh  cells  and  100-200  energy  groups.  These  requirements,  combined 
with  the  probable  necessity  of  running  ~  20  -  40  outer  (energy  group)  iterations  lead  to 
a  highly  computationally  intensive  calculation.  Based  on  these  and  other  investigations 
described  elsewhere  in  this  document,  we  have  concluded  that  the  Monte  Carlo  method 
serves  as  well  or  better  for  calculations  involving  electron-phonon  interactions. 
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m.  Analytical  Benchmark  Verification  of  SNCODE  Algorithm 

1.  Introduction 

The  Sjv  method  is  an  approximate  numerical  method  for  the  solution  of  the 
transport  equation.  Due  to  the  fact  that  we  generally  obtzdned  good  S yv  agreement  with 
standard  Monte  Carlo  calculations  of  energy  deposition  in  solids,  there  was  adequate 
reason  to  believe  that  the  algorithms  employed  in  SNCODE  were  providing  “correct” 
answers.  At  the  least,  this  agreement  provided  indication  that  the  integrals  of  the 
scalar  flux  over  energy  checked  out.  This  agreement  did  not,  however,  reveal  much 
information  pertaining  to  the  accuracy  of  the  flux  calculations.  A  collaborative  effort 
with  B.  D.  Ganapol  and  J.  C.  Garth^^'provided  tin  excellent  opportunity  to  monitor  the 
accuracy  of  the  spatial  differencing  algorithms  of  SNCODE.  Ganapol  developed  a 
multigroup  benchmark  calculation  that  yielded  exact  determination  of  the  group  scalar 
fluxes  for  electrons.  The  availability  of  this  benchmark  method  made  it  possible,  for 
the  first  time,  to  make  absolute  determinations  of  the  accuracy  of  the  multigroup-S;v 
method  as  applied  to  electron  transport.  We  give  here  a  brief  summary  of  the 
methodology  and  report  some  of  the  results  obtained. 

2.  Multigroup  Benchmark  Calculations 

The  multigroup  Boltzmann  equation^^®’^^^  (Eq.  20)  provides  the  theoretical  basis 
for  the  approximate  calculation  of  electron  transport,  scattering  and  energy  loss  in  a 
planar  infinite  medium  (in  practical  terms,  a  medium  larger  than  one  electron  range). 

(''Ie  +  ^  =  , 

£  £  ■  (20) 

o'=i  £=o  'i 

In  the  above,  the  notation  and  symbols  have  their  usual  standard  definitions:  Ax^n)  is 
the  group  g  flux;  is  the  group  mean  free  path;  and  E*  ,  are  the  group  transfer  cross 
sections,  which  are  derived  assuming  applicability  of  contmuous  slowing  down  theory^"l 

As  reported  in  Ref.  9,  Ganapol  derived  an  exact  expression  for  the  Legendre 
moments,  ,  of  the  Fourier  transform  of  the  group  flux  He  then  performs  a 

numerical  inversion  of  the  Fourier  transform  of  the  scalar  group  flux,  <i>g^Q,  to  obtedn  an 
exact  determination  of  the  scalar  flux.  The  transform  of  the  scalar  flux  is  evaluated  by 
solving  a  matrix  equation  involving  sums  over  Legendre  moments.  The  scalar  flux  can 
be  ev^uated  to  arbitrary  accuracy,  depending  on  the  choice  of  the  Legendre  series 
order. 


For  our  benchmark  tests  of  SNCODE,  we  assumed  an  isotropic  source,  located  at 
the  center  of  the  infinite  medium,  of  electrons  in  the  highest  energy  group  {g  =1),  i.e. 


Q(m)  = 


1 

5- 


(21) 


2.i  Results 

The  quantity  taken  as  the  basis  for  comparison  of  the  S;v  benchmark 
calculations  was  the  scalar  flux.  The  approach  taken  in  the  comparisons  between 
multigroup  S ff  and  benchmark  calculations  began  with  a  series  of  one-group  calculations 
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for  the  plane  isotropic  source  embedded  in  an  essentially  infinite  aluminum  medium  in 
which  the  scattering  cross  section  was  taken  to  be  linearly  anisotropic.  For  this  one- 
group  case  we  investigated  the  behavior  of  the  ratio  of  successive  errors  as  the  spaticd 
mesh  size  was  decreased  in  order  to  reproduce  the  known  theory.  The  Syv  and 
benchmark  calculations  were  then  made  for  the  multigroup  cases.  These  calculations, 
carried  out  using  40  uniformly  spaced  electron  energy  groups,  were  made  for  a  200  keV 
electron  source  in  zduminum.  Four  scattering  cases,  corresponding  to  progressively 
higher  scattering  anisotropy  order,  were  investigated.  In  the  first,  an  artificial  isotropic 
scattering  cross  section  was  devised,  and  in  the  remaining  four,  Syv  and  benchmark 
calculations  were  compared  for  extended  transport  corrected  cross  sections  of  order  2 
through  4. 

Because  the  benchmark  calculations  can  only  accommodate  an  infinite  medium, 
and  the  Sjv  algorithm  requires  a  finite  medium,  we  chose  the  dimensions  of  the  test 
mesh  portions  of  the  scattering  media  to  be  deeply  imbedded  within  two  very  thick 
coarsely  meshed  slabs  each  of  thickness  10  mean  free  paths.  For  the  multigroup 
calculations,  we  chose  the  slab  dimensions  to  correspond  to  ±  1  r.u.  (=  range  unit), 
which  for  our  200  keV  electron  source  in  aluminum  translates  into  a  thickness  of 
±0.058  g/cm^.  These  geometries  closely  approximate  infinite  media  for  all  practical 
purposes. 

2.1.1  Plane  Isotropic  Source;  One-Group;  Linearly  Anisotropic  Scattering 

A  set  of  5  Si2  calculations  were  made  with  the  linear  discontinuous  differencing 
algorithm  for  a  plane  isotropic  source  located  at  the  center  (x=0)  of  a  thick  scattering 
medium  extending  20  mean  free  paths  in  both  positive  and  negative  x-directions.  The 
scattering  parameters  were  (  in  units  of  mean  free  path),  Ej  =  1.0,  £«  =  0.8,  = 

E,  /3.  The  spatial  mesh  discretization  step  sizes  for  the  5  calculations  were  taken  to  oe 
=  2*~”;  n=l,2,3,4,5.  In  each  calculation.  Ax  was  held  constant  over  the  spatial 
mesh  test  region  (-10  <x  <  -flO).  Figure  10  is  a  plot  of  the  successive  error  ratios  r„(x) 
[  =  c„(x)/c  ^j(x)].  From  the  curves  shown,  a  region  of  uniformity(3  <  |x|  <  6)  could  be 
identified  tor  successive  discretization  error  ratios.  In  this  region  it  was  found  that  the 
error  always  decreased  for  n=l,2,3.  For  the  »=1  and  n=2  cases,  the  ratio  is  nearly 
constant  and  has  the  value  ~  8,  which  indicates  second  order  accurate  behavior  for  the 
linear  discontinuous  algorithm  (when  the  error  decreeises  as  the  exponent  n,  the 
approximation  is  said  to  be  (n-1)  order  accurate).  Also  in  this  region,  when  n=3,  it  is 
seen  that  the  error  ratio  settles  down  to  ~  8  for  one-half  of  the  interval;  however,  the 
ratio  r4  does  not  show  a  consistent  decrease  in  error  as  Ax  is  decreased  from  0.25  to 
0.125.  This  is  attributable  to  the  fact  that  it  is  not  possible  to  further  improve  the 
result  by  halving  Ax  with  a  fixed  inner  iteration  convergence  precision.  The  fractional 
errors  observed  in  these  cases  were  generally  comparable  to  the  convergence  precision 
(10"®).  Large  variations  in  r„  were  observed  both  near  the  source  plane  and  the  region 
boundary.  The  poor  results  near  the  source  plane  arise  from  the  singular  nature  of  the 
source  geometry,  while  near  the  region  boundary,  changes  in  mesh  size  propagate  errors 
into  the  test  zone.  The  uniform  region  is  situated  well  away  from  both. 
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Fig.  10.  Ratios  of  successive  errors  for  5  one-group  linear 
discontinuous  calculations.  Spatial  disrectization  step  size  is 
varied  according  to  =  n=l-5.  Stable  region,  3<|il<6 

(mean  free  path  units),  is  indicated.^®^ 

2.1.2  Plane  Isotropic  Source;  Multigroup;  Isotropic  Scattering 

Four  sets  of  SNCODE  calculations  were  made  for  the  200  keV  isotropic  electron 
source  imbedded  in  infinite  aluminum.  The  scattering  cross  section  employed  for  these 
calculations  is  physically  unrealistic  in  that  200  keV  electrons  in  aluminum  do  not 
scatter  isotropically.  However  in  this  first  attempt  to  test  the  multigroup  benchmark, 
we  chose  for  convenience  a  cross  section  which  contains  only  the  zero-th  order  Legendre 
coefficient  of  the  P-12  extended  transport  corrected  Mott  cross  section.  These 
calculations  covered  the  following  four  cases:  linear  discontinuous  and  diamond 
differencing,  each  used  with  Sg  and  S12  Gauss- Legendre  quadrature  sets.  In  each  of 
these  four  cases,  5  separate  Syy  calculations  were  made,  corresponding  to  a  graduated 
scale  of  spatial  mesh  resolutions  [Ax  =  .02,  .01,  .004,  .002,  .001  r.u.]  spanning  a  distance 
of  0.2  r.u.  from  the  source  plane  at  x=0.  The  reasons  for  the  restriction  to  0.2  r.u.  are 
twofold:  1)  to  eliminate  vacuum  boundary  effects  in  the  S;v  calculations;  2)  to  operate 
within  the  high  accuracy  range  of  the  analyticcil  benchmark  solution. 

As  will  be  apparent  in  the  graphs  of  Figs.  11  and  12  it  was  found  that  the 
differences  between  Sg  and  Sjj  ejaculations  were  for  practical  purposes  negligible. 
Figure  11  is  a  plot  of  the  scalar  flux  as  calculated  by  the  benchmark  method  (solid 
lines)  and  SNCODE-Sg  (circles)  with  the  diamond  difference  scheme  for  spatial 
resolution  Ax=.004  r.u.  The  10  curves  are  the  energy  spectra  at  x  =  .02(.02).2  r.u.  As 
can  be  seen,  the  agreement  is  excellent.  Figure  12  is  a  plot  for  the  identical  conditions 
except  that  the  S;v  calculation  was  made  with  SNCODE-S12.  While  little  or  no 
sensitivity  to  discrete  ordinate  order  n  was  observed  in  these  two  calculations,  we  did 
find  that  the  accuracy  of  the  SNCODE  calculations  was  highly  sensitive  to  the  spatial 
discretization.  In  Figure  13  we  show  the  same  case  as  given  in  Figure  12  (diamond 
difference,  Si2>  isotropic  scatter),  however  the  spatial  mesh  resolution  was  coarsened  to 
Ax=.01  r.u.  Here  we  observe  significant  departures  from  the  agreement  obtained  with 
Ax=.004  r.u. 
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Fig.  13.  Compajison  of  the  scalar  flux,  Fj(E)  vs.  E 
groups  at  10  x  positions  as  indicated.  Solid  lines 
benchmark  calculation  results;  circles  represent  tl 
diamond  difference  SNCODE-S12  ccilculation. 
interval  Ai  =  0.004  range  units. 
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2.1.3  Relation  Between  Spatial  Resolution  and  Discretization  Error 

To  further  examine  the  sensitivity  of  cailculations  to  spatial  discretization,  we 
made  a  study  of  the  dependence  of  discretization  error  on  resolution.  Having  a  set  of 
true  benchmark  calculations  enabled  this.  Figure  14  is  a  plot  of  the  exponent,  b,  of  the 
relative  error  dependence  [e  oc  (Ax)"^]  of  the  spatial  resolution  vs.  group  midpoint 
energy,  for  6  positions  along  x,  as  the  resolution  is  increased  from  Ax=.02  to  Ax=.01 
r.u.  for  the  linear  discontinuous  algorithm.  The  relative  error  (t)  was  evaluated  as  the 
deviation  of  the  Sjsf  result  from  the  benchmark  result.  The  inner  iteration  convergence 
precision  for  the  upper  graph  was  set  to  10  while  for  the  lower  graph,  we  chose  10”®. 
As  anticipated,  the  exponent  is  -3  as  dictated  by  theory.  It  should  be  noted  that  this 
holds  true  for  a  wider  energy  range  as  the  accuracy  of  the  S;v  calculation  was  increased. 
However,  at  some  point  sufficiently  low  on  the  energy  scale,  no  further  improvement 
was  achieved.  As  might  be  expected,  the  more  accurate  calculation,  to  the  extent 
that  this  crossover  behavior  exists  at  all,  exhibited  this  behavior  at  a  lower  energy 
value.  We  believe  that  this  is  explained  by  an  accumulation  of  the  downscatter  source 
error,  which  is  not  predominant  in  the  high  energy  groups.  It  can  also  be  noted  that 
tightening  of  the  convergence  precision  mitigated  most  of  the  erratic  behavior  of  the 
exponent  at  the  two  positions  (x=.02,  .06  r.u.)  closest  to  the  source  plane. 


Fig.  14.  Plots  of  the  exjwnent  [b]  of  the  relative  error  [t],  where  eaAa:”®,  as  Ai  is 
decreased  from  .02  to  .01  r.u.  The  behavior  of  b  is  shown  as  a  function  of  energy  at  the 
six  X  positions  indicated  for  two  linear  discontinuous  Sg  calculations  with  inner  iteration 
convergence  precisions  given  by  10'^  (upper  graph)  and  10"*(lower  graph)^®^ 
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2.1.4  Plane  Isotropic  Source;  Anisotropic  Scattering 

The  200  keV  electron  source  (at  x=0)  and  aluminum  scatterer  configuration 
(jx|  <  1.0  r.u.)  chosen  for  the  anisotropic  scattering  benchmark  tests  were  the  same  ais 
those  used  in  the  previous  isotropic  scattering  case.  It  had  already  been 
demonstrated^" that  a  realistic  description  of  the  transport,  or  at  least  the  resultant 
energy  deposition,  of  200  keV  electrons  in  aluminum  could  be  achieved  using  a  P-12 
extended  transport  corrected  screened- Rutherford  scattering  cross  section  in  conjunction 
with  the  Sjsf  method.  At  the  risk  of  sacrificing  physical  realism,  we  made  successive 
application  of  the  extended  transport  correction  from  P-1  to  P-4  to  solve  the  “same” 
problem  if  only  approximately.  This  provided  us  with  a  means  for  systematically 
increasing  the  scattering  anisotropy  order  of  the  benchmark  calculations.  Figures  15-17 
show  plots  of  the  scalar  flux  at  10  equispaced  x  values  [0.04(0.04)0.4]  for  the  P-2,  P-3 
and  P-4  extended  transport  corrected  cross  sections  (the  L-th  Legendre  coefficient  of  the 
P-L  extended  transport  cross  section  is  zero,  hence  these  represent  scattering  anisotropy 
orders  1,2,3).  In  all  cases  the  linear  discontinuous  differencing  algorithm  was  employed 
in  40  group  SNCODE-Sjj  calculations.  The  spatial  discretization  step  was  held  constant 
at  Ax=0.002  r.u.  As  before,  the  solid  curves  represent  the  benchmark  solution.  The 
agreement  between  the  two  methods  ranges  between  the  third  and  fifth  significant 
figure  for  all  points  except  at  the  source  plane. 

These  investigations  served  not  only  as  a  demonstration  of  the  first  multigroup 
benchmark  transport  solution,  but  also  as  a  validation  of  the  operation  of  SNCODE.  It 
should  be  pointed  out  that  the  benchmark  algorithm  is  analytically  exact  in  spatial 
dependence  only.  Both  Sj^  and  benchmark  calculations  treat  energy  dependence 
identically.  The  two  calculation  methods,  multigroup- S;v  and  benchmark, 
complemented  one  another.  The  excellent  agreement  between  muItigroup-S/v  and  the 
benchmark  results:  1)  provided  verification  of  the  predicted  error  behavior  for  both  the 
linear  discontinuous  and  diamond  difference  schemes;  and  2)  allowed  us  to  proceed  more 
confidently  to  more  realistic  Syy  electron  transport  cjdculations  with  SNCODE. 
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Fig.  15.  Comparison  of  the  scalar  flux,  Fj;(E)  vs.  E,  for  40  energy 
groups  at  10  x  positions  as  indicated.  Solid  lines  represent  the 
benchmark  calculation  results;  circles  represent  the  results  of  a 
linear  discontinuous  SNCODE-S12  calculation.  Anisotropic 
scattering  of  200  keV  electrons  in  A1  with  P-2  extended  transport 
corrected  cross  sections. 
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Fig.  16.  Comparison  of  the  scalar  flux,  Fj(E)  vs.  E,  for  40  energy 
groups  at  10  x  positions  as  indicated.  Solid  lines  represent  the 
benchmark  calculation  results;  circles  represent  the  results  of  a 
linear  discontinuous  SNCODE-S12  calculation.  Anisotropic 
scattering  of  200  keV  electrons  in  A1  with  P-3  extended  transport 
corrected  cross  sections. 


Aluminum,  200  keV 
P-4  Tronsport  Correction 


10 


.24 


Energy  (MeV) 


Fig.  17.  Comparison  of  the  scalar  flux,  Fi(E)  vs.  E,  for  40  energy 
groups  at  10  x  positions  as  indicated.  Solid  lines  represent  the 
benchmark  calculation  results;  circles  represent  the  results  of  a 
linear  discontinuous  SNCODE  — Sj2  calculation.  Anisotropic 
scattering  of  200  keV  electrons  in  A1  with  P-4  extended  transport 
corrected  cross  sections. 


3.  Analytical  Benchmark  Testing  of  One-Group  SNCODE  Calculations  with  a  Material 
Discontinuity 

Another  series  of  analytical  benchmark  tests  were  conducted  for  the  verification  of 
SNCODE.  The  objective  of  this  research,  a  collaborative  effort  with  B.  D.  Ganapol  and 
J.  C.  Garthl'"],  was  the  development  of  an  analytical  benchmark  calculation  for  one- 
group  neutron  transport  in  an  infinite  medium  with  a  material  discontinuity. 
Participation  in  this  research  effort  provided  an  opportunity  for  using  SNCODE  to 
confirm  the  validity  of  the  analytical  benchmark  solution  and  also  resulted  in  an 
evaluation  of  the  accuracy  of  SNCODE  for  one-group  transport  in  dissimilar  media. 

The  analytical  solution  foi  two  interfacing  half  spaces  of  dissimilar  materia]  was 
obtained  using  a  numerical  Laplace  transform  inversion  which  has  proven  to  be 
successful  for  single  half-space  problems^^'*^.  Isotropic  scattering  w^as  assumed  with  a 
general  spatially  distributed  isotropically  emitting  source  in  the  right  half-space.  The 
solution  developed  in  Ref.  13  was  used  to  benchmark  an  SNCODE  which  was  written 
primarily  for  application  to  electron  transport  where  material  boundaries  have  a 
profound  influence. 

3.1  Uniformly  Distributed  Source 

In  Figure  18  a  comparison  between  the  analytical  benchmark  and  SNCODE  is 
shown  for  S2,  S4,  Sg,  and  S12  calculations.  Isotropic  scattering  was  assumed,  and  the 
albedos  for  the  left-  and  right-hand  half-spaces  were  C2  =  0.8,  Cj  =  0.95,  respectively.  A 
uniform  isotropic  source  is  distributed  throughout  the  right  half-space.  Note  that,  as 
required,  the  flux  for  large  x  approaches  1/(1-  q)  =  20.  For  6  and  12  ordinates,  the 
results  are  virtually  indistinguishable  from  the  analytical  solution.  As  will  be  shown, 
this  close  agreement  did  not  generally  hold  for  sources  located  at  the  interface. 


Fig.  18.  Comparison  of  SNCODE  discrete  ordinates  and 
benchmark  results  (52,54,56  and  5j2)^^^^. 
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3.2  Plane  Isotropic  Source 


Several  test  runs  were  made  for  a  plane  isotropic  source  located  at  the  interface 
between  two  dissimilar  media.  Every  possible  combination  of  Cj  and  for  c  =  .8,  .9, 
.95  and  .99  was  tested.  The  SNCODE  calculations  were  made  for  S2,  S12  and  S2o»  thus 
providing  a  wide  range  of  angular  quadrature  precision.  As  the  order  of  the  S;v 
calculation  was  increased,  the  approximate  result  more  closely  approached  the 
analytical  benchmark,  as  anticipated.  However,  even  the  S20  calculation  did  not  closely 
match  the  benchmark  curve  near  the  material  interface.  Figures  19  and  20  show  the 
plane  source  results  for  C2  =  0.8,  Cj  =  0.95,  and  Cj  =  0.95,  =  0.8,  respectively.  We 

foimd  that,  for  these  calculations,  increasing  the  angular  quadrature  precision  alone 
provided  limited  improvement  in  the  accuracy  of  the  results.  Better  agreement  near 
the  interface  was  obtained  when  the  spatial  discretization  resolution  was  increased  along 
with  the  quadrature  order. 
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Fig.  20.  Comparison  of  SNCODEl-Sj.Sjjj-Sjo  and  benchmark  scalar  flux 
results  as  a  function  of  distance  from  tne  material  interface-source  plane 
(c2  =  0.95,Ci  =0.8)f^^l 

4.  Analytical  Benchmark  Testing  of  One-Group  SNCODE  Calculations  in  Half-Space 
Geometry 

In  addition  to  the  infinite  medium  electron  and  neutral  particle  transport 
investi^tions  described  above,  a  further  collaboration  with  B.  D.  Ganapol  and  J.  C. 
Garth^^mas  recently  been  initiated  and  continues  to  the  present.  This  is  a  series  of 
analytical  benchmark  calculations  were  made  and  then  used  to  test  the  performance  of 
the  Sjv  algorithms  of  SNCODE.  The  analytical  benchmark  calculation  was  formulated 
by  Ganapol  for  one-group  neutral  particle  transport  in  half-space  geometry,  and  is  also 
based  on  the  Laplace  transform  inversion  method^^^l  A  variety  of  source  configurations 
have  been  tested.  Of  these,  the  major  effort  has  been  directed  thus  far  toward 
boundary  source  problems.  An  example  of  our  results  is  given  in  Figure  21.  Here  we 
have  shown  a  benchmark  solution  (solid  curve)  for  the  scalar  flux  in  a  semi-infinite, 
isotropically  scattering  medium  with  albedo  c  =  0.5  where  the  source  is  specified  by  the 
flux  boundary  condition  2  - 1 

=  )  ,  (22) 
where  A  is  the  normalization  constant.  Four  Sjv  cadculations  were  made,  corresponding 
to  n  =  6,  12,  20,  30.  It  is  apparent  from  the  graphs  of  Fig.  21  that  the  and 
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Scalar  Flux 


benchmark  calculations  converge,  as  in  the  previously  presented  results.  This  work  is 
presently  is  progress,  and  benchmark  comparisons  will  be  made  not  only  for  several 
source  configurations,  but  also  for  other  independent  calculation  methods  such  as 
iterative  solutions  of  the  integral  form  of  the  Boltzmann  equation. 
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Fig.  21.  Comparison  of  analytical  benchmark  2ind  S;v  determinations  of 
sc^ar  vs.  I  in  an  isotropically  scattering  medium  (c  =  0.5)  for 

flux  entering  free  surface.  Solid  curve  represents  analytical 
solution;  dashed  curves  represent  SNCODE-56,5j2,52o,-?3o  solutions^^^l 
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IV.  Monte  Carlo  Transport  Calculations  for  Electrons  in  Strong  Electric  Fields  in 
Silicon  Dioxide^'®^ 


1.  Introduction 

Transport  calculations^^®^  for  electrons  in  the  presence  of  a  strong  electric  field  in 
the  oxide  layer  of  an  idealized  MOS  device  were  made  with  the  LOWEND^®^  Monte 
Carlo  code.  In  the  electron  energy  range  8.9  eV  <  E  <20  eV,  where  two  processes, 
electron- phonon  scatter  and  electron-electron  scatter  overlap,  it  was  reasonable  to 
expect  that  a  significant  amount  of  impact  ionization  should  occur,  since  calculated 
electron  energy  spectra  in  the  presence  of  strong  electric  fields  have  high  energy  tails, 
resulting  from  energy  promotion  by  the  electric  field  and  phonon  absorption,  extending 
well  beyond  the  ionization  threshold  (Eg^p  =  8.9  eV).  One  motivation  for  this 
investigation  was  an  attempt  to  explain  the  apparent  inability  to  experimentally 
observe  substantial  secondary  electron  production  in  Si02  and  the  high  energy  tail 
observed  in  the  vacuum  emission  spectr«?^^  ~  .  A  complete  report  of  the  results  of 

these  Monte  Carlo  studies  is  given  in  Ref.  16,  a  collaborative  effort  with  J.  N.  Bradford. 

2.  Discussion 

In  this  work  we  demonstrated  how:  1)  a  revised  calculation  of  the  ionization  cross 
section;  2)  a  change  in  the  value  of  the  deformation  potential  in  acoustic  phonon 
scatter;  and  3)  variation  of  the  electric  field  profile  as  a  result  of  negative  charge 
trapping  in  the  oxide  reduce  the  high  energy  tail  of  the  electron  spectrum.  These  three 
items  combine  to  produce  a  substantial  reduction  in  secondary  electron  production. 
Further,  in  order  to  simulate  experimental  conditions  as  closely  as  possible,  we  added  a 
transport  calculation  through  a  200  A  aluminum  coating  at  the  anode  side  of  the  oxide 
layer.  For  this  we  made  calculations  of  inelastic  electron-electron  scattering  cross 
sections  based  on  the  Tung  -Ritchie^^®^  dielectric  response  function  and  then  wrote  a 
separate  Monte  Carlo  code  incorporating  these  cross  sections. 

A  revision  of  the  impact  ionization  cross  section  curve^^^  for  Si02  was  calculated 
based  on  the  Anderson- Crowell^^^^  ionization  threshold  value  Their  corrected  value 
for  is  11.9  eV,  as  contrasted  with  the  generally  accepted  8.9  eV  value,  and  is 
determined  through  an  cissertion  that  the  ionization  threshold  is  dependent  on  the  ratio 
of  the  electron  effective  mass  values  for  the  valence  and  conduction  bands.  As  can  be 
seen  from  the  two  curves  shown  in  Figure  22,  the  inverse  mean  free  path  (or  cross 
sectio^  for  impact  ionization  is  significantly  lowered  as  a  result  of  the  upwaird  shift  in 
E^f,.  This  combined  with  the  fact  that  ionization  cannot  occur  for  energies  less  than 
11.9  eV  rather  than  8.9  eV,  resulted  in  lowering  the  production  of  secondary  electrons  . 

The  effects  of  a  change  in  the  acoustic  phonon  scatter  deformation  potential  and 
oxide  charging  on  secondary  electron  production  were  studied  with  the  original  impact 
ionization  cross  sections  {E^l^  =  Eg^p).  This  was  done  in  order  to  isolate  these  factors 
from  changes  in  the  ionization  cross  section.  Figure  23  shows  the  fraction  of  secondary 
electrons  produced  at  600A  penetration  depth  in  the  oxide  as  calculated  by  LOWEND. 
The  effects  of  increasing  the  deformation  potential  from  3.5  eV  to  5.0  eV  and  oxide 
charging  are  investigated  separately.  The  charged  oxide  condition  was  simulated  by 
specifying  an  anode  field  of  8  MV/cm  which  decreases  linearly  to  zero  at  the  cathode 
end,  resulting  in  an  average  field  vdue  of  4  MV/cm  in  the  oxide  layer.  Contrasted  with 
this  are  LOWEND  runs  in  which  the  field  value  is  held  const2int  at  4  MV/cm  across  the 
oxide  layer.  This  simulates  an  uncharged  oxide  with  the  same  field  strength  as  the 
charged  oxide  average  field  strength.  The  combined  effect  of  increased  deformation 
potential  and  oxide  charging  is  considerably  lower  secondary  electron  production. 
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Fig.  22.  Change  in  inverse  mean  free  path  for  impact  ionization  according 
to  the  Anderson- Crowell  formula.  Results  for  threshold  energy  of  A/ZE 
as  calculated  by  Ashley^^^’^^l 
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Fig.  23.  Fraction  of  secondary  electrons  produced  by  5000  primary 
electron  case  histories  in  traversing  600A  penetration  depth.  Charged 
oxide  values  are  shown,  by  backward  projection,  to  equal  those  from 
uncharged  oxide  with  field  equal  to  average  of  charged  oxide  field^^®l 
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The  LOWEND  Monte  Carlo  code  was  also  used  to  calculate  the  vacuum  emission 
kinetic  energy  spectra  with  the  5  eV  deformation  potential,  with  and  without  oxide 
charging.  The  spectra  shown  in  Figure  24  were  obtained  with  two  electric  field 
configurations:  linearly  varying  field,  8  MV/cm  at  the  anode  with  average  value  6 
MV/cm  (solid  histogram);  and  constant  electric  field  throughout,  8  MV/cm  (dashed 
histogram).  The  two  cases  shown  illustrate  that  the  value  of  the  field  at  the  anode,  and 
not  the  average  field  value,  determines  the  energy  spectrum  shape.  The  second 
spectrum  graph.  Figure  25,  further  demonstrates  the  important  role  of  the  anode  field 
value.  The  two  spectra  shown  here  result  from  different  anode  field  values,  8  MV/cm 
and  6  MV/cm.  The  average  field  for  both  spectra  is  6  MV/cm.  One  case  (solid 
histogram)  corresponds  to  a  constant  6  MV/cm  field,  and  the  other  (dashed  histogram) 
to  a  linearly  varying  field. 

Finally,  we  performed  a  Monte  Carlo  simulation  of  the  transport  of  electrons 
emitted  from  the  oxide  layer  through  a  200A  thick  aluminum  anode.  The  input 
spectrum  from  the  Si02  was  the  spectrum  shown  for  the  8MV /cm  anode  field  (Fig.  24, 
solid  histogram).  The  inelastic  electron-electron  cross  sections  were  ccilculated  using  the 
electron  gas  model  of  Tung  and  Ritchie^^^^  which  is  outlined  in  Section  2.1. 


Fig.  25.  Comparison  of  kinetic  energy  spectra  with 
same  average  field  values  (6MV/cm)  but  differing 
anode  fields,  6(  — )  and  8( — )  MV/cm^ 
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2.1  Inelastic  E^ectron-Ellectron  Scatter  in  Conducting  Materials 


where 


E  =  kinetic  energy  of  the  incident  charged  particle  (electron), 

Planck’s  constant, 

/ik=  momentum  transfer  (incident  electron  to  target  electron), 
hu  =  energy  transfer  (incident  electron  to  target  electron), 
m  =  incident  electron  rest  mass, 
e  =  electronic  charge, 

e(k,w)  =  exact  dielectric  function  of  the  solid, 
and  k^  correspond  to  the  maximum  (+)  and  minimum  (-)  allowed  momentum 

transfer  vtilues;  ih— 

k±  =  ±  ,  (24) 

For  conducting  materials,  Ritchie^^^  first  proposed  the  use  of  the  Lindhard^^^^form 
for  €(k,w),  the  dielectric  response  function  for  conduction  band  electrons,  as  follows: 

€(k,a;)  =  1  +  (4)  [fi(x,z)  +  if2(x,z)]  ,  (25) 

where  —  e  /irftvp  ,  ^ 

z  =  k/2kF , 

X  =  ftw/  Ep, 

Ep  =  Fermi  ener^  of  the  scattering  material, 

Vp  =  (2  Ep/m)'^  =  velocity  of  electron  with  kinetic  energy  Ep, 

kp=  (2m  Ep/ft^)^'^  =  wave  vector  corresponding  to  the  Fermi  momentum; 


the  functions  fj  and  fj  are  given  by 


fi(x,z)  =  i  +  ^[l-(z-x/4z)2]  +  [1-(z+x/4z)2]  ln||±-^^^i:i||(26a) 

f2(x,z)  =  I 


irx/8z  ; 
ir/8z[  l-(z-x/4z)^]  ; 
0  ; 


Eq.  23  may  be  rewritten  as  2  t  f 

where,  from  Eq.  24,  _ 

z  =  K-Jx  +  1  ±  1)  . 


z  +  x/4z  <  1 
|z-x/4z|<l<z+x/4z 
|z-x/4z|>l 


2  ' 


(26b) 

(27) 

(28) 


In  Eq.  28  we  have  assumed  that  the  incident  charged  particle  is  an  electron.  This 
assumption  will  be  carried  through  the  remainder  of  this  discussion.  Thus,  the  constant 
outside  the  integral  is  written  in  terms  of  the  Bohr  radius  a,,  =  me^/ft^. 


We  evaluated  the  cross  section  function  for  ionization  due  to  excitation  of 
conduction  band  electrons  in  aluminum.  It  was  more  convenient  to  express  the 
scattering  cross  section  not  as  an  explicit  function  of  the  incident  electron  energy,  E, 
and  energy  transfer,  Aw,  but  rather  in  terms  of  the  dimensionless  quantities  f  (not  to  be 
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confused  with  the  dielectric  function)  and  x.  The  parameter  x  has  been  previously 
defined.  The  parameter  e  is  the  incident  electron  energy  (measured  from  the  Fermi 
level)  in  units  of  the  Fermi  energy: 

.  =  (E  -  Ef)/Ep.  (29) 

The  cross  section  for  electron-electron  inelastic  scatter  is  then  given  as 


’■ee(«>x) 


x^Q(f-x)  f  ...  zf^ 

irao(l+f)  J  (z^-|-x2fi)H(y%)^ 

Zi _ 


(30) 


where  ©(t-x)  is  the  step  function  which  insures  that  the  incident  electron  energy 
remains  above  the  Fermi  level. 


A  consequence  of  expressing  the  cross  section  for  electron-electron  scatter  as  in  Eq. 
30  is  that  the  quantity  (1  -}-  t)  r^^(e,x)  is  independent  of  f,  except  for  the  step  function, 
and  could  be  readily  evaluated  by  numerical  integration.  We  have  evaluated  this 
function  for  the  range  0  <  x  <  7.  Our  result,  shown  in  Figure  26,  agrees  with  that 
obtained  by  Tung  and  Ritchie^^*^. 


X 


Fig.  26.  Evaluation  of  normalized  differential  inverse  mean  free 
path[(l -f  f)rgg(f,z)]  for  inelastic  scatter  with  individual  conduction 
band  electrons  in 
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The  cross  sections  shown  in  Figure  26  were  incorporated  into  a  Monte  Carlo  code 
that  we  wrote  for  the  purpose  of  “transporting’  the  electrons  emerging  from  the  Si02 
layer  through  the  aluminum  layer.  The  trajectories  of  the  primary  electrons  and  the 
secondaries  produced  in  the  aluminum  layer  were  traced,  paying  due  attention  to  the 
multiple  boundary  crossings  at  the  oxide-conductor  interface.  The  electrons  that  re¬ 
entered  the  oxide  layer  were  “banked”,  pcissed  through  the  LOWEND  code,  and  re¬ 
transported  in  the  aluminum  as  many  times  as  was  necessary  to  produce  a  stable 
transmission  spectrum  from  the  aluminum.  The  resulting  vacuum  emission  spectrum  is 
shown  in  Figure  27.  The  integral  transmission  fraction  was  0.064.  This  figure  included 
energetic  secondaries  produced  in  the  aluminum.  As  can  be  seen,  there  is  no  high 
energy  tail  as  was  experimentally  observed^’®’ As  stated  in  Ref.  16,  there  is  reason  to 
believe,  based  on  discussions  among  researchers  in  this  field,  that  the  aluminum  layers 
used  in  the  experiments  may  have  contained  pinholes,  and  that  the  experimental  data 
can  be  reproduced  by  linearly  combining  a  proportionately  small  amount 
(fraction  ~  0.08)  of  the  oxide  anode  spectrum  with  the  aluminum  emission  spectrum. 


Fig.  27.  Transmission  spectrum  through  200A  A1  anode  adjoining 
600A  Si02  layer.  Input  spectrum,  Si02-»A1,  is  characterist  ^'  of 
an  8MV/cm  anode  field^^^l 


36 


V.  Monte  Carlo  Calculations  of  the  Transport  of  MeV  Electrons  in  Insulators 
1.  Introduction 


The  recent  availability  of  experimental  measurements  of  electric  fields  and  charge 
distributions  in  such  materials  as  polymethylmethacrylate(PMMA)  by  Zahn,  Wright, 
Hikita,  et  has  afforded  an  opportunity  to  study  charge  build-up  in  dielectric 

materials,  an  important  factor  in  spacecraft  charging.  To  study  this  and  other  data, 
and  to  perform  numerical  simulations  of  dielectric  charging  and  breakdown,  we  wrote  a 
Monte  Carlo  program  ELMCF  to  calculate  current,  charge  and  energy  deposition 
profiles  in  the  presence  of  electric  fields.  Of  existing  electron  Monte  Carlo  codes,  the 
most  widely  used  (and  tested)  is  the  code  series,  which  is  incapable  of  handling 

problems  involving  applied  electric  fields.  We  were  advised  by  one  of  the  ITS  authors 
that  modification  of  this  code  to  incorporate  electric  fields  would  probably  destroy  the 
code’s  utility.  For  this  reason,  we  undertook  to  write  our  own  ELMC,  a  simple  Monte 
Ceirlo  program  that  makes  use  of  the  same  basic  scattering  algorithm  as  is  used  by  the 
ITS  series,  the  multiple  scattering  (condensed  collision)  formulation  of  Berger^^’^’^^l 
After  we  were  satisfied  that  ELMC  performed  correctly,  we  incorporated  electric  field 
effects  into  this  program  md  named  the  expanded  version  ELMCF. 


2.  ELMC  -  Condensed  Collision  Electron  Monte  Carlo  Code 

The  electron  scattering  model  used  in  ELMC,  is  based  on  the  screened-Rutherford 
cross  section.  Straggling  and  secondary  electrons  are  not  considered.  Electron  current, 
energy,  charge  deposition,  transmission  and  reflection  are  calculated  for  a  single 
material  slab  medium.  The  same  algorithms  as  are  found  in  the  ITS  series  are  used 
here  to  cedculate  the  Legendre  expansion  coefficients  of  the  screened-Rutherford  cross 
section  and  the  collision  stopping  power  (including  the  energy-dependent  density 
correction  terms).  The  single  collision  cross  section  Legendre  coefficients  are  then  used 
to  calculate  the  cumulative  Goudsmit-Saunderson  multiple-scattering  angular 
distribution^^^l  The  basic  idea  underlying  the  multiple-scattering  angular  distribution 
is  that  in  the  process  of  slowing  down,  the  electron  traverses  a  known  path  length 
segment,  say  As,  for  a  given  energy  loss  AE.  Along  this  path  length  segment  the 
electron  undergoes  many  single  collisions  with  other  electrons  in  the  scattering  medium. 
Depending  on  the  electron  energy,  a  representative  path  length  segment  may  account 
for  as  many  as  300  to  400  collisions.  For  the  energy  range  of  interest  we  adopted 
Berger’s  aJgorithm^^*^  of  selecting  N  electron  energy  interveils  where  the  upper  bounds 
(E„,  E„^i  .  .  E/v)  of  successive  intervals  (aE„,  AE„^.i,  .  ,  AE/y^),  decreases  by  2”*^®  (= 
.9170).  With  this  as  a  guide,  we  then  computed  the  As„  using  the  continuous  slowing 
down  approximation. 

The  remainder  of  the  Monte  Carlo  program  is  a  simple,  standard  trajectory  analog 
code.  The  number  of  energy  intervals,  N,  (or  pathlength  segments)  depends  on  where 
one  chooses  the  low  energy  cut-off,  E^.,  for  the  calculation.  We  chose  to  use  the  default 
value  used  in  the  ITS  codes  of  E^  =  0.05  Eo,  where  Eg  is  the  source  energy.  This 
condition  translate  into  34  energy  intervals,  when  the  source  energy  Wcis  chosen  as  the 
highest  energy  value  in  the  cross  section  table.  In  the  ELMCF  program,  the  number  of 
energy  intervals  generally  exceeded  34,  since  the  presence  of  an  electric  field  required 
allowance  for  the  possibility  of  energy  promotion  above  the  source  value. 

The  number  of  trajectory  steps  actually  used  in  our  Monte  Carlo  calculations  was 
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Energy  Deposition  (MeV— cm2/g) 


approximately  3  times  the  number  of  energy  intervals.  In  order  to  obtain  an  improved 
representative  sampling  of  the  multiple-scattering  angular  distribution,  we  split  the  As„ 
into  equal  size  sub-steps  within  each  AE„,  in  this  case  3  sub-steps  of  length  As„/3,  and 
then  sampling  the  multiple-scattering  distribution  3  times  as  often.  This  is  the  same 
procedure  as  appccirs  in  the  ITS  codes.  In  the  following  four  figures,  we  show  some  of 
the  energy  and  charge  deposition  profile  results  obtained  with  ELMC,  and  compare 
them  with  results  obtained  with  TIGER,  the  one-dimensional  ITS  Monte  Carlo 
program.  Both  programs  were  run  using  50000  case  histories.  The  first  two  figures 
(Figs.  28  and  29)  compare  the  energy  and  charge  deposition  profiles  in  PMMA  for  the 
case  of  a  0.4  MeV  electron  beam  normally  incident  on  a  slab  of  thickness  0.14  g/cm^, 
which  corresponds  to  ~  1.2  range  units.  The  results  are  not  identical,  but  track  pretty 
well.  The  ITS  codes  contain  a  more  complete  cross  section  library  that  includes 
inelastic  scatter.  We  believe  that  the  differences  are  attributable  to  subtle  differences  m 
the  evaluation  of  the  multiple-scattering  models.  In  the  second  two  figures  (Figs.  30 
and  31)  the  energy  and  charge  deposition  profiles  in  aluminum  are  compared  for  a  5.0 
MeV  electron  beam  incident  on  a  3.2  g/cm^  slab(  ~  1.2  range  units).  As  can  be  seen. 


incident  electron  beam  in  PMMA. 
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Energy  Deposition  (MeV— cm2/g) 
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Fig.  30.  Comparison  of 
deposition  profiles  calculated  by  ITS* 
and  ELMC.  5.0  MeV  normally  incident 
electron  beam  in  Al. 


Fig.  31.  Comparison  of  charge  deposition  profiles 
calculated  by  and  ELMC.  5.0  MeV  normally 

incident  electron  beam  in  Al. 
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Charge  Deposition  and  Current 


In  Figure  32  we  show  current,  energy  and  charge  deposition  profiles  for  one  of  the 
experimental  configurations  reported  in  Ref.  25,  a  2.6  MeV  electron  beam  normally 
incident  on  a  PMMA  slab  of  thickness  1.4  g/cm^.  The  electron  beam  is  first  diffused  by 
0.1  g/cm^  PMMA  before  entering  the  1  A  g/cm^  slab.  The  attenuation  accounts  for  the 
incident  current  having  a  value  slightly  less  than  unity.  The  histogram  curves  of  Fig. 
32  represent  the  first  iteration  (  zero  applied  field)  in  a  series  of  Monte  Carlo 
calculations.  Each  subsequent  iteration  (Monte  Carlo  calculation)  will  be  made  using 
ELMCF  with  internal  electric  fields  calculated  from  the  current  profile  curve  of  its 
predecessor. 
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Energy  Deposition  (MeV— cm2/g) 


We  also  compared  results  obtained  with  ELMC  with  energy  cind  charge  deposition 
profiles  as  calculated  by  our  discrete  ordinates  program,  SNCODE.  In  the  next  two 
figures  (Figs.  33  and  34)  we  compare  the  0.4  MeV  electron  beam  profiles  as  obtained  by 
the  two  methods.  The  physical  models  used  in  both  methods  are  based  on  the  same 
screened-Rutherford  scattering  cross  sections.  ELMC  incorporates  these  into  the 
Goudsmit-Saunderson  multiple-scattering  model,  while  SNCODE  employs  the  P-12 
extended  transport  correction.  The  results  obtained  by  the  two  methods  converge  cis 
the  number  of  energy  groups  (and  the  computational  effort)  is  increased  from  40  to  120. 
Depending  on  the  problem  parameters,  the  ELMC  run  time  on  a  386  PC  for  50000  case 
histories  is  typically  of  the  order  of  1-2  hours. 


Fig.  34.  Comparison  of  charge  deposition  profiles 
calculated  by  ELMC  and  SNCODE(100  and  40  groups). 
0.4  MeV  normally  incident  electron  beam  in  Al. 
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3.  ELMCF-Condensed  Collision  Electron  Monte  Carlo  Code  for  Applied  Electric  Field 
Problems 

Inclusion  of  electric  field  effects  in  our  electron  transport  Monte  Carlo  calculations 
required  modification  of  the  penetration  (PENET)  and  energy  tally  (SCORE) 
subroutines  of  ELMC.  In  the  penetration  routine,  the  principal  modification  was  that  of 
dividing  the  pathlength  increment  into  fine  sub-intervals  5s„(=As„/15)  and  calculating 
the  change  in  electron  trajectory  and  energy  within  that  sub-interval  brought  about  by 
the  acceleration  or  deceleration  due  to  the  presence  of  the  electric  field.  By  the  use  of 
sub-intervals,  it  was  anticipated  that  each  of  these  incremental  changes  would  amount 
to  a  small  perturbation  in  the  overall  trajectory  and  that  their  cumulative  effect  would 
be  the  smooth  incorporation  of  energy  and  trajectory  changes  over  a  path  length  step 
As„.  The  following  two  graphs  (Figs.  35  and  36)  illustrate  the  effects  of  uniform 
positively  directed  (along  z)  applied  electric  fields  on  energy  and  charge  deposition 
profiles.  The  calculations  were  made  for  a  2.6  MeV  electron  beam  normally  incident  on 
PMMA.  There  are  four  histograms  in  each  graph,  corresponding  to  8  =  0.0,  0.5,  1.0 
1.5  MV/cm.  In  the  simulation  of  the  experiment^  situation  described  in  Ref.  25,  the 
maximum  field  strength  encountered  is  ~  1.5  MV/cm  and  is  in  general  not  uniform 
across  the  slab.  ELMCF  has  the  capability  to  operate  with  spatially  varying  electric 
field  8j. 


Fig.  35.  Energy  deposition  profiles,  ais  calculated  by  ELMCF,  in 
PMMA  in  the  presence  of  a  positively  directed  constant  electric 
field  for  a  2.6  MeV  normally  incident  electron  beam.  Results  for  4 
field  values  are  shown  (0,  0.5,  1.0,  1.5  MV/cm). 
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Fig.  36.  Charge  deposition  profiles,  as  calculated  by  ELMCF,  in 
PMMA  in  the  presence  of  a  positively  directed  constant  electric 
field  for  a  2.6  MeV  normally  incident  electron  beam.  Results  for  4 
field  values  are  shown  (0,  0.5,  1.0,  1.5  MV/cm). 
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ROME  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search,  development,  test,  and  technology  transition  in  support  of  Air 

3 

Force  Command,  Control,  Communications  and  Intelligence  (C  I)  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 

O 

ESD  elements  to  perform  effective  acquisition  of  C  l  systems.  In  addition, 
Rome  Laboratory's  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainability  and  testability. 


